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Abstract 



A magnetic bunch compressor was designed and is be- 
ing commissioned to provide higher peak current for the 
Advanced Photon Source's (APS) Low-Energy Undulator 
Test Line (LEUTL) free-electron laser (EEL) Of great 
concern is limiting emittance growth due to coherent syn- 
chrotron radiation (CSR). Tolerances must also be carefully 
evaluated to find stable operating conditions and ensure that 
the system can meet operational goals. Automated match- 
ing and tolerance simulations allowed consideration of nu- 
merous configurations, pinpointing those with reduced er- 
ror sensitivity. Simulations indicate tolerable emittance 
growth up to 600 A peak current, for which the normal- 
ized emittance will increase from 5 to about 6.8 /im. The 
simulations also provide predictions of emittance variation 
with chicane parameters, which we hope to verify experi- 
mentally. 

1 INTRODUCTION 

A companion paper [^] reviews magnetic bunch compres- 
sion and shows a schematic of our system. I assume the 
reader is familiar with this paper The APS bunch com- 
pressor design is an outgrowth of studies [Ij] by P. Emma 
and V. Bharadwaj of Stanford Linear Accelerator Center 
(SLAC). They explored a number of designs, including 
symmetric and asymmetric four-dipole chicanes. Starting 
from this work, I investigated a large number of configu- 
rations with various values of R^q, asymmetry, and final 
current. Eor each configuration, detailed longitudinal and 
transverse matching was performed, followed by tracking 
with CSR and wakefields. Then, sensitivity analysis was 
performed for all configurations, followed by jitter simula- 
tions for the least sensitive configurations. 

This work relied on elegant [^, a 6-D code with a fast 
simulation of CSR effects, plus longitudinal and transverse 
wakefields. elegant also performs optimization of actual 
tracking results, such as bunch length, energy spread, and 
emittance. 

Simulation of the linac uses the RFCA element, a matrix- 
based rf cavity element with exact phase dependence. Our 
linac has quadrupoles around the accelerating structures. 
Hence, I split each 3-m section into about 20 pieces, be- 
tween which are inserted thin-lens, 2nd-order quadrupole 
elements. A series of such elements is used for each 
quadrupole. 

A Green's function technique is used to model wake- 
fields, using a tabulation of the SLAC-structure wake func- 
tions provided by P. Emma |pl]. To reduce running time. 



one longitudinal wake element is used per 3-m section, 
which is a good approximation for relativistic particles. For 
transverse wakes, I used one wake element per rf cavity el- 
ement (about 20 per section). 

The CSR model used by elegant is based on an equa- 
tion for the energy change of an arbitrary line charge 
distribution as a function of the position in the bunch and 
in a bending magnet. Details of this model will be pre- 
sented by the author at an upcoming conference. Effects of 
changes in the longitudinal distribution within a dipole are 
included. CSR in drift spaces is included by propagating 
the terminal CSR "wake" in each bend through the drifts 
with the beam. 

2 MATCHING 

Longitudinal and transverse matching has the goal of pro- 
viding configurations for the 300-A and 600-A LEUTL op- 
erating points [^. The starting point for the simulations is 
macro particle data generated [0| with PARMELA, giving 
the 6-D distribution after the photoinjector (PI). See Figure 
1 in [|]. 

Longitudinal matching involves adjusting the phase and 
voltage of L2 (see [^] for nomenclature) to obtain the de- 
sired current and energy after the chicane. Then, L4 and 
L5 are adjusted to minimize the energy spread and obtain 
the desired final energy. Longitudinal matching includes 
longitudinal wakefields, rf curvature, and higher-order ef- 
fects in the beam transport, by matching tracked properties 
of the simulated beam. 

Figure |l] shows the longitudinal phase space for the 300- 
A case with R^q = —65 mm, which exhibits a current spike 
of nearly 1200 A. The matching ignores this spike (which 
is shorter than a slippage length for 530 nm) because of the 
way "current" is defined [||]. 
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Figure 1 : Typical longitudinal phase space (300-A case) 




Figure 2: Typical twiss parameters in the chicane region. 



Following longitudinal matching, transverse matching is 
done for each configuration. Initial Twiss parameters are 
obtained from the rms properties of the PARMELA beam. 
Starting values for the quadrupoles were obtained from 
matching "by hand" for one configuration. Four sequen- 
tial elegant runs work the beta functions down the linac. 
The most important constraints maintain small beta func- 
tions in the linac (for transverse wakefield control), small 
horizontal beta in dipole B4 (to reduce CSR effects), and 
matching for the emittance measurement sections. Figure 
^ shows sample Twiss parameters in the chicane region. 

The matching is highly automated, so that only the de- 
sired beam current and energy needs to be specified. Eval- 
uation of tolerances and randomized simulations are also 
automated, being set up by scripts from the corresponding 
matching runs. Transfer of data between simulation stages 
is handled using SDDS files and scripts ||^], reducing errors 
and increasing the number of configurations that can be ex- 
amined. For example, a script is used to scan all configu- 
rations and give power supply specifications. A distributed 
queue utilizing 50 workstations is used to run the simula- 
tions. 

Figure |] shows emittance vs. R^e for the symmetric 
(A=l) and asymmetric (A=2) cases at 300 A and 600 A. 
For 300 A, the symmetric and asymmetric cases are very 
similar. For 600 A, the difference is 10% or more, which 
should be measurable. 

One surprise in Figure || is that the emittance does 
not uniformly increase as |i?56| increases, even though 
elegant shows the expected monotonic increase (due to 
CSR) vs bending angle for a single dipole with a constant 
input beam distribution. This is apparently due to variation 
in the compressed bunch distribution between cases with 
the same "current" but different R^q. For smaller |i?56|, 
there are higher current spikes at the head of the bunch, 
leading to a larger and more rapidly changing CSR wake, 
which in turn leads to larger emittance growth. The effect 
is even more pronounced in the 1200- A cases (not shown). 
Insertion of the scraper between B2 and B3 to remove the 
low-energy part of the beam can reduce the height and 
width of the current spike, resulting in lower emittance. 
Unfortunately, this also reduces the current in the rest of 



the bunch considerably. 

Earlier simulations showed that emittance trends can be 
changed significantly by inconsistent values of the hori- 
zontal beta function at the exit of B4. All of these sub- 
tleties will make for difficult interpretation of experiments 
in which R^q is varied. However, because compression to 
different currents for fixed i?56 involves only adjustment 
of the rf phases and voltages, comparision of the emit- 
tance growth for different amounts of compression should 
be more straightforward. 
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Figure 3: Horizontal normalized emittance vs. Rqq. 



3 TOLERANCE DETERMINATION 

Tolerances are driven by the FEL gain length, trajectory, 
and wavelength stability requirements [^. The 10% rms 
gain length variation limit is easy to use in elegant as it 
computes FEL performance directly using Xie's parame- 
terization [p^. Beam trajectory limits (^50 /im, ^50/ir) 
are included separately as they are not incorporated into 
Xie's formula. The 1-nm rms wavelength variation limit 
is a challenging goal at 530 nm as it puts a 0.1% limit on 
energy variation. 

The analysis begins by running single-parameter 
"sweeps" to assess the effect on the constrained quanti- 
ties (gain length, trajectory, and wavelength) of accelera- 
tor parameters (e.g., rf phase). Sweeps included rf phase 
and voltage; photoinjector timing, charge, and energy; and 
chicane dipole strength. From these sweeps, a script deter- 
mines the limit on each parameter change due to the var- 
ious specifications, showing that configurations with the 
largest i?56 are least sensitive to difficult-to-control tim- 
ing and phase errors. These configurations experience the 
most emittance degradation from CSR, but tend to yield the 
shortest gain length as they have the smallest energy spread 
(L2 being closer to crest). 

The limits, shown in Table |l], are larger for the 600-A 
case because the 1-nm wavelength constraint is easier at 
120 nm than 530 nm. Nine parameters are limited primar- 
ily by the wavelength constraint and four others by hori- 
zontal trajectory constraints. Hence, to determine the rms 
tolerance, one simply divides each sweep limit by VlV, N 



Table 1: Selected sweep limits for R^q = — 65 mm 
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17° 
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L4/L5 phase 


0.77° 


1.45° 


L2 voltage 


0.11% 


0.31% 


L4/L5 voltage 


0.52% 


1.4% 


PI timing 


0.29 ps 


0.88 ps 


PI energy 


0.26% 


1.1% 


PI charge 


12% 


>20% 



Table 3: Results of 300 randomized simulations with re- 
laxed tolerancelevelsfori^ 





300 A 


600 A 


quantity 
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% 
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% 




jitter 
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jitter 
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(a;) (iim) 


89 


72 


89 


81 


{xl) {fir) 


59 


64 


68 


58 


(y) {^rn) 


63 


88 


127 


79 


{y) (m^) 


138 


62 


245 


39 


Lgain (m) 


0.048 


68 


3 


1.3 


X (mn) 


9.6 


9 


2.8 


27 



being the number of parameters limited by a particular con- 
straint. For the horizontal trajectory, iV was doubled to 
eight to allocate half the budget to nonswept parameters 
(e.g., corrector magnets). Some of these phase and timing 
tolerances are beyond the state of the art. 



4 RANDOMIZED SIMULATIONS 

Randomized simulations were used to confirm the toler- 
ances and examine errors not covered by the sweeps (e.g., 
corrector jitter, quadrupole jitter, and alignment). These 
were done for the most stable configurations (i.e., R^q = 
—65 mm). Because some tolerances are beyond the state 
of the art, I used randomized simulations to determine the 
impact of "relaxed" tolerances, assuming these rms levels 
[pl|: 1° rf phase jitter, 0.1% rf voltage jitter, 1 ps timing 
jitter, 5% charge jitter, and 2% PI energy jitter 

Tables || and |^ show the results, respectively, for the 
sweep-derived tolerance levels and the relaxed levels. The 
sweep-derived tolerance levels result in meeting the speci- 
fications for the FEL, while the relaxed levels, not surpris- 
ingly, do not. One surprise in the relaxed case is the large 
jitter in the vertical plane. This results from uncorrected 
nonlinear dispersion in a vertical dogleg between the linac 
and the LEUTL, a problem which can be readily remedi- 
ated using two sextupoles [Bt]. 



Table 2: Results of 300 randomized simulations with 
sweep-determined tolerance levels for R^q = —65 mm 
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0.83 
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